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The electron self-energy of the ti u states in A3C60 (A=K, Rb) is calculated using the so-called 
GW approximation. The calculation is performed within a model which considers the t\ u charge 
carrier plasmon at 0.5 eV and takes into account scattering of the electrons within the t\ u band. A 
moderate reduction (35 %) of the tin band width is obtained. 



I. INTRODUCTION 

The alkali doped fullerenes A3C60 (A=K, Rb) have a 
low-energy charge carrier plasmon at an energy of about 
0.5 eV,ErH which essentially results from the oscillation of 
the alkali electrons donated to the t\ u band. This plas- 
mon has a rather strong coupling to the electrons and it 
is believed to play an important role for the anomolously 
large band width seen in angular integrated photoemis- 
sion for these systems ,□ since it can cause large binding 
energy plasmon satellites. 

It is interesting to calculate the electron self-energy 
due to the coupling to the plasmons, to more generally 
study the effects of the plasmons on the electronic prop- 
erties. We here use the so-called GW approximation, 
where the self-energy is calculated to lowest order in the 
screened Coulomb interaction.El This approximation has 
been widely applied to the electron gas, free-electron like 
system and semiconductors. For these systems it is found 
that the self-energy has a moderate effect on the elec- 
tron dispersion and that it can somewhat increase or re- 
duce the effective mass, depending on the system. A3C60, 
however, has a quite different band structure from these 
systems, and it is interesting to ask if the effective mass 
and the width of the band might be more strongly mod- 
ified in this system. 

This issue is, for instance, raised by the results for the 
optical conductivity. Experimentally, the Drude peak 
is very narrowjj perhaps an order of magnitude nar- 
rower than predicted by a one-particle theoryfl This may 
suggest that the effective mass is greatly enhanced in 
A 3 C6o- On the other_hand, the specific heat has been 
estimated to be small,u perhaps even smaller than would 
have been expected from band structure calculations,u 
suggesting that the effective mass is not strongly renor- 
malized, and perhaps even somewhat reduced. The cou- 
pling to phonons with the coupling constant X p h reduces 
the dispersion by a factor of 1/(1 + X p h)- If this result is 
directly taken over for the plasmons, we may estimate a 
coupling constant A p / ~ 2.5 and a substantial reduction 
(more than a factor of three) of the dispersion. 

The coupling to a plasmon studied here corresponds 
to the GW method in the plasmon pole approximation, 
which was first introduced and extensively studied by 



Lundqvist et a/.Ei and later by Overhausei-0 for the elec- 
tron gas. Shirley and Leuie have applied the GW approx- 
imation to solid C6o-li3 Since they considered undoped 
Ceo, the t\ u plasmon considered here did not enter their 
calculation, while the contribution to the self-energy in 
their calculation (a plasmon at about 25 eV) is not in- 
cluded in the model used below. The two calculations 
are therefore complementary and address different parts 
of the physics of A^Cqq. 

The coupling to the plasmon leads to a so-called plas- 
maron satellite in the_spectral function, where a plas- 
mon has been excited, ill This satellite is obtained in the 
GW approximation)!!] but a more accurate treatment 
would aise-give higher satellites with several plasmons 
excited.EJ'li-3 These higher plasmon satellites are believed 
to be important for the broad spectimn seen in angular 
integrated photoemission for A 3 C6o-ra!j The GW approx- 
imation can therefore not be expected to give an accu- 
rate spectrum in the satellite region. We also observe 
that the Coulomb interaction U between two electrons 
on the same Cgo pplecule is large compared with the 
t\ u band width W)¥i which suggests that many-body ef- 
fects beyond the GW approximation may be important. 
For instance, a proper treatment of the three-fold degen- 
eracy of the t\ u orbital is important for understanding 
why these systems are not Mott-JJubbard insulators, in 
spite of the large value of U/Wx3 The GW results for 
A3C60 should therefore be treated with a certain caution, 
but it is still interesting to see if a strong renormalization 
of the band width is obtained. 

In Sec. II we present the formalism and in Sec. Ill 
the model used in these calculations. The results are 
presented in Sec. IV and discussed in Sec. V. 



II. FORMALISM 

We write the dielectric function of A 3 Ceo as 
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where eo is the contribution to e due to all excitations 
in the system except the ones inside the t\ u band. Since 
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most of these excitations have a rather high energy rel- 
ative to the energy scale of interest here (the t\ u band 
width), we assume eo to be energy independent. The sec- 
ond term describes the excitations inside the t\ u band. 
u)Qp is the plasmon frequency one would deduce for free 
electrons with the same density as the t\ u electrons and 
m* is the band mass of the t\ u electrons. The value of us 
where e(q, us) = gives the plasmon frequency 
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Above we have neglected local-field effects as well as 
the q-degendence of eo- Neither approximation is quite 
justifiedE Thus the q-dependence of e tends to give lu p 
a positive dispersion, while the local-field effects tend to 
give a negative dispersion. As far as the plasmon fre- 
quency is concerned, however, these two effects essen- 
tially cancel, in agreement with the experimental obser- 
vation that us p has a negligible dispersion.cl Eqs. (0J!) 
therefore give a good description of the plasmon fre- 
quency, but neglect the substantial broadening of the 
plasmon.B 

For K3C60 the t\ u electron density corresponds to the 
electron gas density parameter r s = 7.3ao and usq p = 2.4 
eV. In an electron gas of this density, the occupied part of 
the band is about 0.9 eV, while thje calculated full band 
width of K3C60 is about 0.6 eVfl This corresponds to 
to* ~ 0.9/(0.6/2) = 3. At q = £ ~ 4.4E3 We then 
find that us p ~ 0.66 eV. This is somewhat larger than the 
experimental result, oj p ~ 0.5 eV,0 which will be used in 
the following. 

■The self-energy is in the GW approximation written 
asE 
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while the second part gives a correlation contribution. 
We can interpret 2us p /(us 2 — cj|) in Eq. (f|) as a boson 
Green's function, and 
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as a coupling constant, in analogy with previous work.0 
Closing the integration contour in the upper half of the 
complex us plane, we obtain the correlation contribution 
to the self-energy. 
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III. MODEL 



We consider a model with three t lu orbitals. The hop- 
ping matrix elements between these orbitals. hav .e been 
described in a tight-binding parametrization,r 9 tr°l which 
is used here. This parametrization is used to calculate 
the noninteracting band structure energies £kn and wave 
functions 
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where 
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where U nn » n > n » (k — k ) is the Coulomb matrix element 
between the Bloch states |kn) and |k n ) with the argu- 
ment r and |kn ) and |k n ) with the argument r . Here 
k is a wave vector and n a band index. + is an positive 
infinitesimal number, which assures the proper behavior 
for large us . £ k '„" is the noninteracting energy of the 
state |k n ), and fio is an average of the self-energy over 
the Fermi surface. We split the dynamically screened 
Coulomb interaction U/e in two parts 
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is a Bloch state of a t% u molecular orbital & u (r). There 
are N molecules with the coordinates Rj 

We further have to specify the matrix elements of the 
Coulomb interaction. We assume 
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For the Coulomb matrix elements between Bloch states 
we then find 



The first term then gives a statically screened exchange 
contribution to the self-energy. 
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We replace the sum over R by an integral over all space 
outside a Wigner-Seitz sphere with the radius Rq = 5.56 
A. This gives 

e 2 
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where f2 is the volume of the Wigner-Seitz cell. We have 
put Uq = 4 eV, using a simple estimate based on the 
radius of the Cgo molecule. Alternatively, we can extend 
the integration over R over all space, putting R = 0. In 
this case we should put Uq = to avoid double counting. 
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where £kn is the band structure energy. We can also 
obtain the quasi-particle weight as 

1 
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where the derivative is evaluated at the quasi-particle en- 
ergy. The resulting quasi-particle (GW) band structure 
is compared with the noninteracting tight-binding (TB) 
and Hartree-Fock (HF) band structures in Fig. El. The 
figure shows a moderate reduction (~ 35%) of the GW 
band width relative to the TB width, while the HF band 
width is almost a factor of two larger (75%) than the TB 
result. 

The quasi-particle strength is relatively small, Z ~ 
0.4 — 0.5. This is smaller than what is found in electron 
gas calculations at metallic densities. Since the three t\ u 
electrons correspond to a very low density with r s ~ 7, 
this small value of Z is, however, not very surprising. Ac- 
tually extrapolation of earlieia GW calculations for the 
electron gas to r s = 7.3 suggest Z ~ 0.52. The small 
value of Z implies that a substantial spectral weight is 
transferred to satellites. 

Replacing Uq — 4 eV and Rq = 5.56 A in the model 
( |l2| ) by Uq — and Rq = 0, e.g., a pure 1/q 2 interaction 
gives a small change of the quasi-particle band structure. 
This suggests that the results are not very sensitve to the 
details of the model. 



V. DISCUSSION 
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FIG. 1. Quasi-particle (GW) (a) band structure for the 

t lu band compared with the noninteracting tight-binding 
(TB) (b) and Hartree-Fock (HF) (c) band structures. Rela- 
tive to the TB band structure, the GW band width is reduced 
by about 35 %, while the HF width is increased by about 75 
%. 



The results in Fig. [| illustrate that there is a large can- 
cellation between exchange and correlation effect. This 
is also observed in the electron gas, although-in-that case 
the cancellation is more complete than hereoEll 

Of particular interest is the energy dependence of the 
self-energy. For the sake of the discussion, we neglect the 
q-dependence of the coupling constant g(q) (Eq. (^)) and 
replace it by its average over the Brillouin zone. We fur- 
ther assume that the plasmons couple to a nondegenerate 
band of width 2B and with a constant density of states. 
The correlation part of the self-energy is then given by 
(Eqs. ©) 
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IV. RESULTS 
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We have calculated the self-energy as described above. 
The quasi-particle energy E^n is obtained by solving the 
equation 



is a coupling constant defined as for the electron-phonon 
coupling and N(Q) = 1/(2B) is the density of states. 
In the present case we have A ~ 2.5. In the limit 
| a; | << Lu p << B we have 
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In the opposite limit where uj p >> B, \u>\ we have 
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The coupling constant (g/cup 



1.4 in the present case. 



The "phonon"-like formula in Eq. (17]) predicts that the 
quasi-particle weight is Z ~ 0.29 while the formula ( |lS| ) 
predicts Z ~ 0.42. The latter result fits the actual calcu- 
lations nicely, not too surprisingly in view of the parame- 
ter range considered here (ui p = 0.5 eV, B = 0.3 eV). The 
fact that the system is not in the "phonon" parameter 
range (oj p << B) therefore means that correlation effects 
lead to a less drastic narrowing of the band. This nar- 
rowing is furthermore to a substantial part compensated 
by exchange effects. 

Nevertheless, the narrowing of the band is larger than 
found in the GW approximation for the electron gas at 
metallic densities. We want to discuss this further. 
We rewrite the self-energy as 
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In the first term the sum is over the whole Brillouin zone, 
and it can therefore be shown that within the model in 
Eq. ( |l2| ) it is state independent. We can then focus on 
the next two terms. If u> p is very large, these two terms go 
to zero. The band is then just shifted to lower energies, 
without any change in shape or width. We next consider 
finite values values of lo p , but for a moment we assume 
that Up is still larger than the band width. If u + /iq is 
put at the bottom of the band, both terms are then posi- 
tive for all k . States at the bottom of the band are then 
pushed upwards by the last two terms. In the same way, 
states at the top of the band are pushed downwards, lead- 
ing to a narrowing of the band. These arguments are not 
qualitatively different if u p is somewhat smaller than the 
noninteracting band width. Some of the energy denomi- 
nators in Eq. (19) could become negative, but the effect 
of reducing oj p tends to be a further reduction of the band 
width. We observe that these argument are specific for 
this model, where we have assumed that the states have 
not just a lower bound but also an upper bound. For the 
electron gas and most other systems there is no upper 
bound, and no definite statements of this type can be 
made. 

In the calculation above, we have only considered the 
coupling to the t\ u charge carrier plasmon and neglected 



the coupling to, e.g., the plasmon at about 25 eV as well 
as the exchange interaction with all occupied states ex- 
cept the t\ u states. These effects -mere considered in the 
calculation by Shirley and Louie,E3 who found a broad- 
ening of the t\ u band by about 30 %. If this broadenig 
is added to our results we find an essentially unchanged 
t\ u band width. This result can then not explain the 
narrow Drude peak in the optical conductivity, but it is 
essentially consistent with the rather small specific heat 
deduced for these systems. 
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